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Abstract 

We consider non-supersymmetric two-dimensional CP(N — 1) mo¬ 
del deformed by a term presenting the bosonic part of the twisted 
mass deformation of J\f = 2 supersymmetric version of the model. Our 
deformation has a special form preserving a Z]\r symmetry at the La- 
grangian level. In the large mass limit the model is weakly coupled. Its 
dynamics is described by the Higgs phase, with Z jv spontaneously bro¬ 
ken. At small masses it is in the strong coupling Coulomb/confining 
phase. The Z jv symmetry is restored. Two phases are separated by 
a phase transition. We find the phase transition point in the large-A 
limit. It lies at strong coupling. As was expected, the phase transition 
is related to broken versus unbroken Zn symmetry in these two re¬ 
spective phases. The vacuum energies for these phases are determined 
too. 


1 Introduction 


As well known, two-dimensional CP(N — 1) model is an excellent theoretical 
laboratory for modeling, in a simplified environment, a variety of interesting 
phenomena typical of non-Abelian gauge theories in four dimensions [1, 2], 
Recently, two-dimensional CP(N — 1) model was shown to emerge [3] as 
a moduli theory on the world-sheet of non-Abelian flux tubes presenting 
solitons in certain four-dimensional Yang-Mills theories at weak coupling 
[4, 5, 6, 7]. The flux tube solutions in the bulk (“microscopic”) Yang-Mills 
theory depend on an adjustable parameter of dimension of mass. When this 
parameter is large the flux tubes are in fact Z N strings; they evolve towards 
non-Abelian strings as the above mass parameter decreases and eventually 
vanishes. Correspondingly, the world-sheet theory is not just the CP(N — 1) 
model; rather it is the CP(N — 1) model mass-deformed in a special way 
that preserves a Z N symmetry of the model. The mass term we deal with 
coincides with a special choice of the twisted mass [8] in supersymmetric 
CP(N — 1) model. We hasten to emphasize that we will focus exclusively on 
non-supersymmtric version, to which the two-phase phenomenon we study is 
inherent. There is no such phenomenon in supersymmetric version. 

In the limit of vanishing mass deformation, the CP(N— 1) model is known 
to be a strongly coupled asymptotically free held theory [9]. A dynamical 
scale A is generated as a result of dimensional transmutation. However, at 
large N it can be solved by virtue of l/N expansion [1], The solution found by 
Witten exhibits a composite photon, coupled to N quanta n, each with charge 
1/VN with respect to this photon. In two dimensions the corresponding 
potential is long-range. It causes linear confinement, so that only n*n pairs 
show up in the spectrum. This is the reason why we refer to this phase as 
“Coulomb/confining.” In the Coulomb/confining phase the vacuum is unique 
and the Zn symmetry is unbroken. 

On the other hand, if the mass deformation parameter is 3> A, the model 
is at weak coupling, the held n develops a vacuum expectation value (VEV), 
there are N physically equivalent vacua, in each of which the ZN symmetry 
is spontaneously broken. We will refer to this regime as the Higgs phase, 
although this name has a Pickwick sense. Usually the Higgs mechanism 
implies that a gauge boson eats a would-be Goldstone meson thus acquiring 
a mass that screens long range interactions. In our case, at m 3> A, there is no 
gauge boson to begin with (see Sect. 4). However, the long-range interaction 
inherent to the Coulomb/confining phase does not take place; that’s why it 
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is not unreasonable to refer to the phase as the Higgs phase. 

In Ref. [3] it was argued that the twisted mass deformed CP(N — 1) 
model undergoes a phase transition when the value of the mass parameter is 
~ A. The argument was largely based on analysis of the flux tubes and their 
evolution in the underlying four-dimensional theory. In this paper we will 
show, basing our consideration on two-dimensional model per se, in the large 
N limit, that a phase transition between Z N broken and unbroken (i.e. the 
Higgs and Coulomb/confining) phases does indeed occur at m 2 = A 2 . The 
change of regimes takes place in a narrow interval m 2 — A 2 = 0(1/N) where 
the method we use is insufficient to resolve details of the phase transition. 
In particular, the task of finding a conformal held theory emerging at the 
critical point m 2 = A 2 remains open. 

The issue of two phases and phase transitions in related models was pre¬ 
viously addressed by Ferrari [11, 12]. While the first paper [11] deals with 
CP(N — 1) models, neither the methods used nor results have a significant 
overlap with the results reported below. In [12] Ferrari exploits 1/N ex¬ 
pansion methods which are similar to ours. The model to which Ref. [12] 
is devoted is a mass-deformed O(N) model with a Z 2 symmetry at the La- 
grangian level. The point of the phase transition separates the Z 2 broken 
phase at weak coupling from the the Z 2 unbroken phase at strong coupling, 
which is in parallel with our result. The phase transition in [12] is argued to 
be of the Ising-model type. We do not expect this to be valid in the Z N case 
we deal with in the present paper. 

The paper is organized as follows. In Sect. 2 we introduce the twisted mass 
deformation of the non-supersymmetric CP(N — 1) model that preserves Zn 
at the Lagrangian level. We also review some well-known facts regarding this 
model in the strong and weak coupling regimes. The weak coupling regime 
(large twisted masses) is especially simple since here the theory can be treated 
perturbatively and exhibits a Higgs-like behavior. At strong coupling we are 
guided by Witten’s large-W solution. In Sect. 3 large-W methods are used 
to solve the model at arbitrary m. The critical point separating Zjy broken 
and unbroken phases is determined and the vacuum energies are calculated 
for both phases. Section 4 presents a remark concerning inadequacy of the 
large-iV methods for determination of the nature of the critical behavior. 
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2 Twisted-mass deformed CP(N — 1) model 
and its phases 

In this section we describe a twisted mass deformation of the CP(N — 1) 
model preserving Z at. Then we discuss its two distinct phases in two opposite 
limits, m ■C A and m 3> A. 

As was mentioned, the origin of the word “twisted” lies in supersymmetry, 
more exactly, extended Af = 2 supersymmetry. Aspects of the supersymmet¬ 
ric version were analyzed by Dorey [10], who found an exact solution in the 
holomorphic sector. We study the non-supersymmetric version, obtained by 
discarding the fermion sector. 

2.1 The model 

As usual in two dimensions, the Lagrangian can be cast in many different 
(but equivalent) forms. For our purposes the most convenient formulation is 
in terms of the n fields. 1 To set our notation, let us first omit the twisted 
mass. Then the CP(N — 1) model can be written as 

S — J d 2 x {( d a + iA a ) n* e ( d a — iA a ) n e + A (n* e n £ — r)} , (1) 

where rr is an IV-component complex hied, t = 1,2subject to the 
constraint 

n* e n e = r (2) 

where r is the inverse coupling constant of the model. More exactly, the 
standard relation between r and g 2 is 

r = 2/g 2 . 

The action (1) and other similar expressions below are given in the Euclidean 
space. The constraint (2) is implemented by the Lagrange multiplier A in 
Eq. (1). The held A a in the Lagrangian is auxiliary too; it enters with no 
derivatives and can be eliminated by virtue of the equation of motion, 

A a = d a n l . (3) 

_ 2 r 

1 They are referred to as “quarks” or solitons in Ref. [1]. 
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Substituting Eq. (3) in the Lagrangian, we rewrite the action in the form 


S — I d 2 x < d a n* e d a ri' + —(n* e d a n e ) 2 + A ( 


n*«n e — r ) 


) • 


( 4 ) 


The model (4) is a generalization of the 0(3) sigma model. The latter is 
formulated as 

r 


S = - / d 2 xd a Sd a S 


(5) 


where S' is a three-component vector subject to the constraint S 2 = 1. At 
N = 2 the CP(N — 1) model (i.e. CP{ 1)) reduces to 0(3) through the 
substitution 

S a = -(n*T a n), a = 1,2,3, (6) 

r 

where r a are the Pauli matrices. The constraint S 2 = 1 follows from Eq. (2), 
while 


d a Sd a S 


1 


£\2 


d a n* e d a n + ~(n* e d a n ) 


(7) 


The coupling constant r is asymptotically free [9], and defines the dy¬ 
namical scale of the theory A through 


A 2 = M 2 v exp 


47rr c 

’Iv" 


( 8 ) 


where M uv is the ultraviolet cut-off and tq is the bare coupling. The combi¬ 
nation N/r is nothing but the ’t Hooft constant which does not scale with 
N. As a result, A scales as N° at large N. One can also introduce the 9 
term, if one so desires, 

Se = ^— j d 2 x£ ai d a A 1 = J d 2 xe ai ( d a n* t <9 7 n £ ) . (9) 

Now let us add to the action (1) or (4) a mass term of a special form, 

S rn = f d 2 x ~ me)™*} , (10) 

J i 

where o is an auxiliary complex held (with no kinetic term), and we choose 

f 2 j7ti a, \ 

rrii — m exp f J , ^ = 0,1,..., N — 1. (11) 
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The parameter m in Eq. (11) can be assumed to be real and positive. This 
is not the most general choice of the twisted mass deformation. In general, a 
single condition is imposed, = 0; which destroys SU(1V)/U(1) preserv¬ 

ing only residual U(l) Ar_1 . We want to maintain an additional Zn symmetry, 
however, which is automatic under (11). The Zjy symmetry of the action has 
historic roots [3], but what is important at present is that the twisted mass 
deformed model with Z^ symmetry is interesting on its own. The Z jv sym¬ 
metry plays an important role in identifying a phase transition between the 
Higgs and Coulomb/confining phases of the theory. 

Eliminating cr by virtue of the equation of motion, 

a = - mi n*i n e , (12) 

r i 

we get 

f 

72 2 J 1 2 ^ u * t 1 

a x m < r -> e N n» n > . 

\ r r j 

It is instructive to see what becomes of this mass term at N 
0(3) sigma model). Then, Eq. (13) implies 

S m — r J d 2 xm 2 (1 — Sf), (14) 

which is obviously Z 2 symmetric. If m is large, m S> A, the theory has two 
vacua, at S :i = 1 and S :i — — 1. In both vacua there are two elementary 
excitations, Si and S 2 , with masses 2m. 

In the general N case the action at hand has the following Zn symmetry: 

• 27rfc 

cr —> e n a, n e -> n t+k 

for every fixed l and k — 1,2,..., N. (15) 

2.2 The Higgs phase 

At large m, m S> A, the renormalization group flow of the coupling constant 
is frozen at the scale m. Thus, the model at hand is at weak coupling 
and the quasiclassical analysis is applicable. The potential (10) have N 



(13) 

= 2 (i.e. the 
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degenerate vacua which are labeled by the order parameter (a), the vacuum 
configuration being 

a = rri ( 0 , n e ° = \fr , and n e = 0 if £ ^ £ 0 . (16) 

In each given vacuum the Z^ symmetry (15) is spontaneously broken. 

There are 2(IV — 1) elementary excitations 2 with physical masses 

M e = | m e - m io |, £^£ 0 . (17) 

Besides, there are kinks (domain “walls” which are particles in two dimen¬ 
sions) interpolating between these vacua. Their masses scale as 

Mf nk ~rM £ . (18) 

The kinks are much heavier than elementary excitations at weak coupling. 
Note that they have nothing to do with Witten’s n solitons [1] identified as 
solitons at strong coupling. The point of phase transition separates these 
two classes of solitons. 

2.3 The Coulomb/confining phase 

Now let us discuss the Coulomb/confining phase of the theory occurring at 
small m. As was mentioned, at rri — 0 the CP(N — 1) model was solved 
by Witten in the large- N limit [1]. The model at small m is very similar 
to Witten’s solution. (In fact, in the large-A" limit it is just the same.) In 
Sect. 3.3 we present a generalization of Witten’s analysis which we will use 
to study the phase transition between the Zn asymmetric and symmetric 
phases. Here we just briefly summarize Witten’s results for the massless 
model. 

If m = 0, classically the field rr can have arbitrary direction; therefore, 
one might naively expect spontaneous breaking of S\J(N) and the occurrence 
of massless Goldstone modes. Well, this cannot happen in two dimensions. 
Quantum effects restore the full symmetry making the vacuum unique. More¬ 
over, the condition (2) gets in effect relaxed. Due to strong coupling we have 
more degrees of freedom than in the original Lagrangian, namely all N fields 
n become dynamical and acquire masses A. 

2 Here we count real degrees of freedom. The action (1) contains N complex fields r/. 
The phase of r/° can be eliminated from the very beginning. The condition njrd = r 
eliminates one more field. 
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k=0 k=l k=0 


Figure 1: Linear confinement of the n-n* pair. The solid straight line rep¬ 
resents the ground state. The dashed line shows the vacuum energy density 
(normalizing £ 0 to zero). 


This is not the end of the story, however. In addition, one gets another 
composite degree of freedom. The U(l) gauge field A a acquires a standard 
kinetic term at one-loop level, 3 of the form 

NA~ 2 F a pF a g. (19) 

Comparing Eq. (19) with (1) we see that the charge of the n fields with respect 
to this photon is 1 /a JN. The Coulomb potential between two charges in two 
dimensions is linear in separation between these charges. The linear potential 
scales as 

V(R) ~Ar (20) 

where R is separation. The force is attractive for pairs n and n, leading to 
the formation of weakly coupled bound states (weak coupling is the manifes¬ 
tation of the 1/N suppression of the confining potential). Charged states are 
eliminated from the spectrum. This is the reason why the n fields were called 
“quarks” by Witten. The spectrum of the theory consists of nn- “mesons.” 
The picture of confinement of n’s is shown in Fig. 1. 

The validity of the above consideration rests on large N. If N is not large 
the solution [1] ceases to be applicable. It remains valid in the qualitative 
sense, however. Indeed, at N — 2 the model was solved exactly [13, 14] 
(see also [15]). Zamolodchikovs found that the spectrum of the 0(3) model 
consists of a triplet of degenerate states (with mass ~ A). At IV = 2 the 
action (4) is built of doublets. In this sense one can say that Zamolodchikovs’ 
solution exhibits confinement of doublets. This is in qualitative accord with 
the large -N solution [1], 

Inside the nn mesons, we have a constant electric field, see Fig. 1. There¬ 
fore the spatial interval between n and n has a higher energy density than 
the domains outside the meson. 

’’By loops here we mean perturbative expansion in 1/IV perturbation theory. 
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Modern understanding of the vacuum structure of the massless CP(N— 1) 
model [16] (see also [17]) allows one to reinterpret confining dynamics of 
the n fields in different terms [18, 3]. Indeed, at large N, along with the 
unique ground state, the model has ~ N quasi-stable local minima, quasi¬ 
vacua, which become absolutely stable at N — oo. The relative splittings 
between the values of the energy density in the adjacent minima is of the order 
of 1/N, while the probability of the false vacuum decay is proportional to 
N^ 1 exp(-N) [16, 17]. The n quanta (n quarks-solitons) interpolate between 
the adjacent minima. 

The existence of a large family of quasi-vacua can be inferred from the 
study of the 6 evolution of the theory. Consider the topological susceptibility, 
i.e. the correlation function of two topological densities 

Jd 2 x(Q(x), Q( 0)), (21) 


where 

Q=± £aid «£i = -j- e ai (cPn* c>V) . (22) 

The correlation function (21) is proportional to the second derivative of the 
vacuum energy with respect to the 6 angle. From (22) it is not difficult to 
deduce that this correlation function scales as 1/N in the large N limit. The 
vacuum energy by itself scales as N. Thus, we conclude that, in fact, the 
vacuum energy should be a function of 9/N. 

On the other hand, on general grounds, the vacuum energy must be 
a 27r-periodic function of 6. These two requirements are seemingly self- 
contradictory. A way out reconciling the above facts is as follows. Assume 
that we have a family of quasi-vacua with energies 


E k {9) ~ NA 2 


1 + const 


/ 2nk + 9 

V 



k = Q...,N -1 (23) 


A schematic picture of these vacua is given in Fig. 2. All these minima 
are entangled in the 6 evolution. If we vary 6 continuously from 0 to 2n 
the depths of the minima “breathe.” At 9 = n two vacua become degenerate, 
while for larger values of 9 the former global minimum becomes local while the 
adjacent local minimum becomes global. It is obvious that for the neibohring 
vacua which are not too far from the global minimum 


Ek+i ~ E k 
8 


TV' 


rs j 


(24) 




Figure 2: The vacuum structure of CP(N — 1) model at 9 = 0. 


This is also the confining force acting between n and n. 

One conld introduce order parameters that would distinguish between 
distinct vacua from the vacuum family. An obvious choice is the expectation 
value of the topological charge. The kinks n e interpolate, say, between the 
global minimum and the first local one on the right-hand side. Then n’s 
interpolate between the first local minimum and the global one. Note that 
the vacuum energy splitting is an effect suppressed by 1/N. At the same 
time, kinks have masses which scale as N °, 

M n kink _ A _ ( 25 ) 

The multiplicity of such kinks is N [19], they form an IV-plet of SU(iV). This 
is in full accord with the fact that the large- N solution of (1) exhibits N 
quanta of the complex held n. 

In summary, the CP(N — 1) model in the Coulomb/confining phase, at 
small m, has a vacuum family with a fine structure. For each given 9 (except 
6 = 7T, 37T, etc.) the true ground state is unique, but there is a large number 
of “almost” degenerate ground states. The splitting is of the order of A 2 /N. 
The Zjy symmetry is unbroken. The spectrum of physically observable states 
consists of kink-anti-kink mesons which form the adjoint representation of 

SU (IV). 

At large m the theory is in the Higgs phase; it has N strictly degenerate 
vacua; the Z^ symmetry is broken. We have N — 1 elementary excitations 
rv with masses given by Eq. (17). 

Thus we conclude that these two regimes should be separated by a phase 
transition [3]. This phase transition is associated with the Zn symmetry 
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breaking: in the Higgs phase the Z^ symmetry is spontaneously broken, 
while in the Coulomb phase it is restored. For N = 2 we deal with Z 2 which 
makes the situation akin to the Ising model. 

3 Solution at m ^ 0 in the large -N limit 

In this section we will generalize Witten’s analysis [1] to include m ^ 0. The 
twisted mass deformed action is 

S— /d 2 J|V^| 2 + A(|nT-r o)+ E|(.-m^| 2 }, (26) 

were V„ = d a — iA a and rrq is defined in Eq. (11), and r 0 is the bare coupling 
constant. 

3.1 Effective theory 

As soon as the action (26) is quadratic in n e we can integrate over these fields 
and then minimize the resulting effective action with respect to other fields. 
Large -N limit ensures that corrections to the saddle point approximation are 
small. In fact, this procedure boils down to calculating one-loop graphs with 
holds n ( propagating inside loops. 

In the Higgs phase the held n e ° develop a VEV. One can always choose 
£ 0 = 0 and denote rr° = n. The held n, along with er, are our order parame¬ 
ters that distinguish between the Coulomb/conhning and the Higgs phases, 
see (16). 

Therefore, we do not want to integrate over n a priori. Instead, we will 
stick to the following strategy: we integrate over N — 1 helds n e with £ ^ 0. 
The resulting effective action is to be considered as a functional of n, A and 
a. To find the vacuum configuration, we will minimize the effective action 
with respect to n, A and a. 

Integration over rr with l ^ 0 produces the determinant 

N—\ 

JJ [det (-<9 2 + A + \a - m e | 2 )] _1 , (27) 

i =i 

where we dropped the gauge held A a . In principle, as was explained in 
Sect. 2.3, quasi-vacua in the Coulomb/co nhning phase have non-vanishing 
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expectation values of the operator (22). However, we cannot see these VEVs 
in the leading order in N. Since the analysis we carry out applies to the 
leading order in the large-TV limit, we set A a = 0. 

Calculating (27) we get the following contribution to the effective action: 


N -1 


47T 


Y (A+ \a-m e \ 2 ) 


t=\ 


M 2 

ln A + |< r - V m ,|2 + 1 


(28) 


where we dropped a quadratically divergent contribution which does not 
depend on A and a. 

Equation (8) implies that the bare coupling constant ro in (26) can be 
parameterized as 


N . 

= ~r ln 

47r 


M 2 

_ UV 

A 2 


(29) 


Substituting this expression in (26) and adding (28) we see that the term 
proportional to A ln M 2 V is canceled out, and the effective action is expressed 
in terms of the renormalized coupling constant, 


T ren 


1 

47T 


N -1 




A + |cr - me 
A 2 


(30) 


where we neglect 0(1/N) contributions. In addition to the coupling constant 
renormalization we have to carry out renormalization of the held a leading 
to a renormalization of its vacuum expectation value. To this end we add 
the corresponding counterterm to the bare action (26), namely, 



where c is a finite constant to be fixed below. This counterterm ensures 
that the infinite term proportional to |cr — me\ 2 ln M 2 V in the determi¬ 

nant (28) is canceled and, the renormalized VEV of a is finite. We fix the 
coefficient c below demanding that (a) —m 0 = 0 in the Higgs phase, see (16). 

Assembling all contributions together we get the effective action in the 
form 


S = 



n 


2 


11 



TV—1 


— ^(A+|TT-m £ | 2 ) 


e =i 

TV-1 


1 — In 


A + | cr — 
A2 


121 


iiEi"- m 'i : 


( 32 ) 


^=1 


Now, minimizing this action with respect A, n and o we arrive at the following 
set of equations: 


— ^ren ? 

(A + |cr - : m 0 | 2 ) n = 0, 


(33) 

(34) 


N—l 


— ^ (ff - m t ) In 


£=1 


A + | cr — rri(_ | 5 
A2 


IV 


+ (a - m 0 ) |n| 2 + — ca = 0 , 


47T 


(35) 


where r ren is given in Eq. (30), and we take into account the fact that the sum 
1 m? is relatively suppressed: instead of 0(N ) it is 0(1) so that we lose 
the factor of N. Equations (33), (34) and (35) represent our master set that 
determines the vacua of the theory. Note that Eq. (33) is a renormalized 
version of the bare condition (2). In addition to the above equations the 
vacuum configuration must satisfy two extra constraints, 


r 


> 0 , 


(36) 


and 

Re A > 0. (37) 

The Erst condition just follows from |n| 2 > 0, see (33). The second one 
becomes clear if we examine the original bare action of the model. From 
Eq. (26) we conclude that the integral over A runs along the imaginary axis 
(remember, we use the Euclidean formulation of the theory.) The saddle 
point solution for A can have (and will have) a non-vanishing real part. To 
ensure convergence of the path integral over rr’s the real part of A at the 
saddle point should be non-negative. It is important that a is an independent 
integration variable, and the integral over rr’s must be convergent at all 
values of a, in particular at a = me. 
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We will see shortly that the constraints (36) and (37) are important con¬ 
ditions which single out physical phases existing in the given range of the 
parameter m. In the subsequent sections we will study solutions to our 
master-set equations and show that at large m the theory is in the Higgs 
phase while at small m it is in the Coulomb/confining phase, the boundary 
being at m — A. 


3.2 The Higgs phase at large N 

At m > A the solution to Eqs. (33), (34) and (35) has the following form: 

(A) = 0, 

(a) = m 0 , 

(n) = y/nZ, (38) 


where we use the gauge freedom of the original model to choose n real, as was 
explained in Sect. 2.2. We see that the fields a and n have non-vanishing 
VEV’s and, as a result, the Z^ symmetry is spontaneously broken. Our 
choice n = no was of course arbitrary. In fact, we have N strictly degenerate 
vacua as shown in Eq. (16). Equation (35) must be used to fix the value of 
the constant c, 


1 




| m e - m 0 
A 2 


(39) 


Substituting this value in the effective action (32), together with VEV’s (38), 
we get the vacuum energy in the Higgs phase, 


-^Higgs vac 0 m . (40) 

2 7T 

The logarithmic term in the second line cancels the third line. 

Now, let us have a closer look at the additional constraints (36) and (37). 
The latter condition is trivially satisfied while to examine the impact of the 
condition (36) we substitute cr = mo = m and A = 0 in the expression (30) 
for the renormalized coupling constant. Then we get 


T rpn — 


— V'hi ^ ~ 

47T ^ 

"=i 


A 2 


1 2m sin ^ 

— > In 

7 T ' 


7r 


i=l 


A 


N m 

= — In — 

2vr A 


(41) 
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where the sum over £ is calculated in the large N limit. The constraint (36) 
implies 

m > A . (42) 

The Higgs phase has a clear-cut meaning at large m. Hence, the above result 
is compatible with intuition. We will see momentarily that the lower bound 
of the allowed domain, m = A, is the phase transition point. 

3.3 The Coulomb/confining phase at large N 

At small m the appropriate solution of the master equations (33), (34) and 
(35) has the form 

o = 0 , 
n — 0 , 

A = A 2 — m 2 . (43) 

The vacuum expectation value of the n held vanishes, as one would expect 
from the Zw symmetric phase, while Eq. (33) is satisfied because 

Tren = 0 (44) 

in the vacuum (43), cf. Eq. (30). In fact Eq. (43) becomes an m ^ 0 
generalization of Witten’s saddle point condition which was used to determine 
VEV of A in [1], Upon consulting with Eq. (26) we conclude that in our saddle 
point the mass of the rv quanta is A, independent of the value of the mass 
deformation parameter m. Indeed, the mass squared —> A + |cr — m# | 2 = A 2 . 
Although this statement might seem counter-intuitive, it is correct. We will 
comment on that in the end of this section. Since both a and n do not 
condense in this Coulomb/conhning vacuum the Zjy symmetry is unbroken. 
The bare condition (2) is relaxed due to (44). The solution exhibits more 
degrees of freedom than are present in the Lagrangian. 

Let us turn now to constraints (36) and (37). The Erst one is satisfied 
trivially while the second one implies 

m < A. (45) 

We see that at m < A the theory is in the Coulomb/co nhning vacuum 
(43) while at m > A it is in the Higgs vacuum (38). The value 

m* = A (46) 
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Figure 3: Normalized vacuum energies (4:TtE VSI jC /N A 2 ) versus m 2 / A 2 . The 
solid line shows the actual vacuum energy, while dashed lines correspond 
to a formal extrapolation of the Higgs and Coulomb/confinement vacuum 
energies to unphysical values of m below and above the phase transition 
point, respectively. 


is the phase transition, or critical point. 

Let us calculate the vacuum energy in the Coulomb phase. Substituting 
the vacuum values (43) in the action (32) and using expression (39) for the 
value of the constant c we arrive at the vacuum energy 

N f Tfl? 1 

^Coulomb vac = — < A 2 + + HR 2 log > , (47) 

where the sums over t are calculated in the large -N limit. 

We plot the vacuum energies (47) and (40) for the Coulomb/conhning 
and Higgs phases as a function of m 2 in Fig. 3. At the point of the phase 
transition (46) energy densities of both phases coincide. Moreover, their first 
derivatives with respect to m 2 at this point coincide too. The Higgs curve, 
naively extrapolated below the phase transition, runs below the Coulomb 
curve which might lead one to conclude that the system always stays in the 
Higgs phase. However, the conditions (36) and (37) produce constrains (42), 
(45) which tell us that at m < A the system is in the Coulomb/co nhning 
phase while at m > A it is in the Higgs phase. 

One can check our results for the vacuum energies in both phases per¬ 
forming the calculations in a slightly different form, through the trace of the 
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energy-momentum tensor. The vacuum energy can be obtained as 


E-iran 



( 48 ) 


where 




N—l 

M uv d Muv + Y, (rriidm t + m}d m *) 
1=0 


£(M UV , mg) , 


(49) 


and £(M UV , mi) is the ultraviolet-regulated Lagrangian of the model. Taking 
into account the classical contribution and the quantum anomaly, we pre¬ 
cisely reproduce the vacuum energies in the Coulomb/confining and Higgs 
phases quoted above. 

To reiterate, at large m, at weak coupling, we have N strictly degenerate 
vacua; the Zn symmetry is broken. At small m, at strong coupling, a mixing 
between these vacua takes over, and N vacua split (see Sect. 2.3). The order 
parameter which marks these vacua is the VEV of the operator (22) which is 
non-zero for exited “vacua” with k ^ 0. In the leading order in N to which 
we are limited, we do not see this vacuum splitting. Our result exhibits a 
single vacuum (43). Moreover, we cannot say anything as to the nature of the 
phase transition. The answer can be found upon inspection of a narrow strip 
\m 2 — A 2 1 ~ 0(1/N) (see Sect. 4) which would require tools going beyond 
those exploited here. 

ft is curious to note that the 6 dependence of physical quantities, albeit 
suppressed at large N, is suppressed differently above and below the critical 
point. If in the Higgs phase the suppression is expected to be exponential, it 
is power-like in the Coulomb/confining phase. 

Finally, we would like to comment on the independence of the n-quanta 
mass on m in the Conlomb/conhning phase. The crucial observation is that 
in the absence of n VEVs the twisted mass term (13) is actually quartic in the 
n fields, rather than quadratic. Therefore, while it contributes to interactions 
of the n fields, its contribution to the n mass cannot appear at order N°. 


4 Can we describe critical behavior? 

The full solution of the phase transition problem requires establishing a con¬ 
formal held theory which governs dynamics at the critical point, ft is no 
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accident that so far nothing has been said regarding this issue. In this sec¬ 
tion we argue that the proper description of the critical behavior would re¬ 
quire methods going beyond the 1/N expansion on which we rely. Thus, this 
question remains open. 

To understand the nature of the phase transition one must identify states 
that become massless at the critical point m = A. Let us undertake this 
endeavor, approaching the critical point from the Higgs side. 

In the Higgs phase, at large m, the theory is weakly coupled, and all 
excitations are massive. There are N — 1 complex degrees of freedom. The 
lowest singularity in the correlation function (A a (x), Hg(0)) is a two-particle 
cut, so that no stable held playing the role of a photon exist. On the other 
hand, at m = 0, there are N quanta of the U(l)-charged n fields, and a 
massless photon. It is natural to ask whether a light photon emerges at 
strong coupling as one approaches the critical point from above. 

To address this issue we modify our effective action (32) including the 
gauge held with a kinetic term induced at one loop as in [1], 

Sffiggs = J d 2 x l-^~ F ad + | V Q n| 2 + ^ (|n| 2 - r ren ) 2 + E mggs vac | , 

where r ren is given in Eq. (41). The kinetic term for the gauge held is induced 
through a loop of the quanta. This loop converges in the infrared domain 
and, as we will see shortly, is saturated by the lightest n e quanta, i.e. £ ~ 1. 
In addition to this kinetic term, we also include a quartic term for n (we 
remind that n = n°) which comes from integration over A in (32) around the 
corresponding saddle point at A = 0 (in the quadratic approximation). The 
corresponding “coupling constants” denoted by e 2 en and e 2 en , (in fact, they 
are momentum dependent; hence, the quotation marks) are 



127T _ 2 4n 

6ren = W)' 


(51) 


where 


N -1 

m = £ / 

o _i J 


[k 2 + |me - m 0 \ 2 ][(k - p ) 2 + \me - m 0 \ 


and p is the external momentum, p a <-> id a . 

The lightest states in the sum (52) correspond to l ~ 1. Their masses 
scale as 

9 

Tfi 

\me - m 0 1 2 ~ vry at £ ~ 1. (53) 
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As long as the gauge field A a and the scalar n are much heavier, they cannot 
be treated as stable point-like bound states. 

As we reduce m, the gauge field A a and the scalar n become lighter 
and eventually may cross the threshold and become genuinely stable bound 
states. To evaluate their masses let us take the low-energy limit in (52), 
assuming that p 2 is much less then the masses of lightest elementary states, 
see Eq. (53). Keeping only £(0) we get from (50) the masses of the gauge 
field (the massive gauge field has one real degree of freedom) and the scalar 
n (the phase of n is eaten by the Higgs mechanism), respectively, 


ren ^ r en| p 2 = o ) 


E(0) = 0.17 


N 2 
2 m 2 ’ 


m n — 2 Ken C ren p2=0 ) 


(54) 

(55) 


where we calculated the sum over £ in (52) at p = 0 numerically in the limit 
of large N. The renormalized coupling constant r ren is given in Eq. (41). 

As m approaches m* = A from above, the coupling constant r ren tends to 
zero and, seemingly, so do the masses of the gauge field A a and scalar n, 


rn„ 


m. r 


P m ' 


(56) 


where 8m — rn — m*. However, these conclusions would be correct only if 
the masses of these bound states were much smaller than the masses of the 
lightest rr quanta. Comparing with (53) we see that this would require 


8m 1 
— -C —. 
A N 


(57) 


In other words, the gauge field A a and the scalar n could become light only 
in a very close vicinity of the critical point. Unfortunately, in the domain 
(57) the expansion in 1/At, on which we heavily rely, explodes, and we cannot 
trust our analysis. 

Summarizing, in the narrow strip (57) near the critical point where light 
composite states could occur - those which could become massless at criti¬ 
cality - the 1/At expansion fails. As a result, we cannot derive the conformal 
field theory which would describe our system at criticality. 
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5 Conclusions 


The mass deformed non-supersymmetric two-dimensional CP(N — 1) model, 
with a special Zn preserving twisted mass term, surfaced recently in connec¬ 
tion with non-Abelian strings in four-dimensional gauge theories [3]. This 
model turns out to be very interesting on its own, as a theory with two 
distinct phases and a critical point at strong coupling. Using the large-A" 
expansion we confirmed the fact of the phase transition in m, determined 
the position of the critical point and calculated the vacuum energies in the 
Higgs and Coulomb/confining phases. The major unsolved problem is deter¬ 
mination of the conformal field theory governing dynamics of the model at 
criticality. 

The use of the large -N expansion allowed us to bypass such question as 
“what particular aspect of the strong coupling dynamics is responsible for 
the change of regimes at m = m* = A?” Although no definite answer to 
this question can be given at the moment, it is tempting to conjecture that 
the phase transition is clue to the fact that at m < m* melted instantons 
of Fateev et al. [20] play a crucial role while at m > m* instantons are 
“individualized” and suppressed. 

Instantons are not suppressed at large N at rn — 0 due to a large en¬ 
tropy factor. The theory can be rewritten as a massive fermion theory [20] 
or, equivalently, as the affine Toda theory at fixed coupling constant. Hence 
it is natural to ask how our large -N one-loop calculation captures nontrivial 
instanton effects. A possible answer can be inferred from the supersymmetric 
version. Supersymmetric theory can be treated in two different ways [21], 
Within the first approach one exploits a similar one-loop calculation, while 
within the second approach summation over nonperturbative configurations 
yields a twisted superpotential of the affine Toda type. The mirror symme¬ 
try of supersymmetric version is responsible for equivalence of the vacuum 
structure in both approaches. In our non-supersymmetric version there is no 
evident notion of the mirror symmetry. However, one could still hope that 
the relation between instanton calculus and the one-loop calculation in the 
linear gauged formulation of CP(N — 1) works in a similar manner. 

The phase transition at some value of the twisted mass we demonstrate in 
the present paper in CP(N— 1) seems to be a more general phenomennon tak¬ 
ing place in a class of asymptotically free non-supersymmetric sigma models. 
In particular such phase transition could be expected in the Grassmannian 
sigma models as well as for many toric target manifolds. To an extent, these 
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phase transitions can be considered as non-supersymmetric counterparts of 
the curves of marginal stability in supersymmetric versions. 
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